Phonon-mediated sticking of electrons at dielectric surfaces 
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We study phonon-mediated temporary trapping of an electron in polarization-induced external 
surface states (image states) of a dielectric surface. Our approach is based on a quantum-kinetic 
equation for the occupancy of the image states. It allows us to distinguish between prompt and 
kinetic sticking. Because the depth of the image potential is much larger than the Debye energy 
multi-phonon processes are important. Taking two-phonon processes into account in cases where 
one-phonon processes yield a vanishing transition probability, as it is applicable, for instance, to 
graphite, we analyze the adsorption scenario as a function of potential depth and surface temperature 
and calculate prompt and kinetic sticking coefficients. We find rather small sticking coefficients, at 
most of the order of 10 -3 , and a significant suppression of the kinetic sticking coefficient due to a 
relaxation bottleneck inhibiting thermalization of the electron with the surface at short timescales. 

PACS numbers: 52.40.Hf, 73.20.-r, 68.43.Mn 



I. INTRODUCTION 

A complete kinetic modeling of atmospheric^, inter- 
stellar— or man-made bounded plasmas^— requires 
boundary conditions for the distribution functions of the 
relevant plasma species (electrons, ions, neutrals), that 
is, a quantitative microscopic understanding of the ele- 
mentary processes at the plasma boundary. Of particular 
importance is the build-up of a quasi-stationary nega- 
tive surface charge, which not only depletes the electron 
density in front of the boundary (sheath formation) but 
also acts as an electron reservoir for surface-supported 
electron-ion recombination and secondary electron emis- 
sion which in turn affect the charge balance in the bulk 
of the plasma^- Despite its unquestioned importance, 
little is quantitatively known about the microphysics of 
electrons at plasma boundaries. It is only until recently 
that we proposed that the charging of plasma boundaries 
can be perhaps microscopically understood in terms of an 
electronic physisorption process i 14 ' 15 

The physisorption scenario applies to a plasma elec- 
tron approaching a metallic or a dielectric boundary pro- 
vided its kinetic energy is large enough to overcome the 
Coulomb barrier due to the charges already residing on 
the boundary but small enough to make the surface ap- 
pear as having a negative electron affinity If the electron 
can convert its energy into internal energy of the bound- 
ary, via exciting elementary excitations of the solid, it 
may get stuck (adsorbed) at the boundary. Later it may 
desorb again if it gains enough energy from the boundary. 

Like physisorption of neutral particles^— physisorp- 
tion of electrons is the polarization-induced temporary 
binding to a surface. It can be characterized by a des- 
orption time and a sticking coefficient. At first glance 
physisorption of electrons seems to be not much differ- 
ent from physisorption of neutral particles. There are 
however important qualitative differences which warrant 
a separate theoretical investigation. 

First, albeit not in the focus of our investigation, the 
long-range 1 / z-tail of the image potential leads to a finite 



electron sticking coefficient at vanishing electron energy 
and surface temperaturei 26 ' 31 This is in contrast to the 
quantum reflection, that is, the vanishing sticking coef- 
ficient, one finds in this limit for the short-ranged sur- 
face potentials typical for physisorption of neutral parti- 
clesM 

Second, the surface potential in which physisorption of 
electrons occurs, in particular at plasma boundaries, con- 
sists of a polarization-induced attractive part and a repul- 
sive Coulomb part due to electrons already adsorbed on 
the surface. The limit of vanishing coverage, very often 
adopted in the theoretical description of physisorption of 
neutral particles, is thus only applicable to the very first 
(last) electron approaching (leaving) the boundary. 

Third, in contrast to physisorption of neutral parti- 
cles, physisorption of electrons has to be always described 
quantum mechanically because the image potential varies 
on a scale comparable to the thermal de Broglie wave 
length of the electron^ This is also the case for ph- 
ysisorption of positronium^r— 

Finally, and this will be the theme of our investiga- 
tion, the polarization-induced image potential supports 
deep states, in addition to shallow ones. Direct transi- 
tions from the continuum to deep bound states are very 
unlikely. Hence, a modeling in terms of a quantum- 
kinetic rate equation for the occupancy of the bound sur- 
face states ; 21 ' 22 and Brenig's distinction between prompt 
and kinetic sticking^ is vital for a correct description of 
electron physisorption. For phonon-controlled adsorption 
and desorption, as it occurs at dielectric surfaces, deep 
states also imply that multi-phonon processes have to 
be taken into account in the calculation of state-to-state 
transition probabilities. This can be done either via an 
expansion of the energy dependent T-matri x 16 ! 18 ' 25 , the 
method we are using^, or via a Magnus-type resumma- 
tion of the time-dependent scattering operator^r— 

In the following we investigate adsorption of an elec- 
tron to a dielectric surface at finite temperature assuming 
an acoustic longitudinal bulk phonon controlling electron 
energy relaxation at the surface. To avoid complications 
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due to finite coverage we focus on the first approaching 
electron. Using the quantum-kinetic rate equation for 
the occupancy of the image states of our previous work 35 
(thereafter referred to as I), where we studied desorption 
of an image-bound electron from a dielectric surface, we 
calculate prompt and kinetic sticking coefficients. Com- 
pared to semiclassical estimates^ they turn out to be 
extremely small. Instead of the order of 10 _1 we find 
them to be at most of the order of 10 -3 . We also analyze 
in detail the adsorption scenario as a function of surface 
temperature and potential depth. Most notable, our re- 
sults reveal an energy relaxation bottleneck prohibiting, 
on a short timescale, thermalization of the electron with 
the surface, that is, the trickling through of the electron 
from upper to deep states. The reduced accessibility of 
deep states makes the kinetic sticking coefficient much 
smaller than the prompt sticking coefficient in contrast 
to what is usually found in physisorption of neutral par- 
ticles.^ 

The remaining paper is structured as follows. In Sec. 
II, we specify the quantum-kinetic approach of our pre- 
vious work concerning desorption (paper I2 5 -) to the sit- 
uation of adsorption and introduce prompt and kinetic 
sticking coefficients. We then briefly recall in Sec. Ill 
the calculation of the state-to-state transition probabili- 
ties based on a microscopic model for the electron-surface 
interaction and an expansion of the T-matrix for the dy- 
namic part of that interaction. Mathematical details not 
given can be found in I. Finally, in Sec. IV, we present 
and discuss our results before we conclude in Sec. V. 



II. ELECTRON KINETICS 



continuum to bound surface states. They give the prob- 
ability for initial trapping, which, according to Iche and 
Nozieres^ and Brenig~ is the first out of three stages of 
physisorption. The second stage encompasses relaxation 
of the bound state occupancy and the third stage is the 
desorption of the temporarily bound particle. 

The second stage, which also includes transitions back 
to the continuum, cannot be captured by simple state- 
to-state transition probabilities. Instead, a quantum- 
kinetic rate equation for the time-dependent occupancy 
of the bound surface states n n (t) has to be employed . 18 i 21 
This equation describes processes on a timescale much 
longer than the lifetime of the individual surface states 
but shorter than the desorption time t 18 ' 21 ' 22 Following 
Gortel and coworkers , 18 i 22 

^n n (t) =^T / [W nn >n n >(i) - W n ' n n n {i)] 



J2 W knU n {t) + J2 T t W nk jk{t) , (3) 



where W n 'n is the probability for a transition from a 
bound state n to another bound state n', Wkn an d W nk 
are the probabilities, respectively, for a transition from 
a bound state n to a continuum state k and vice versa, 
and 



jk(t) = n k (t)r t 1 



(4) 



is the incident electron flux which in principle can be 
non-stationary. 

The solution to Eq. ([3]) can be obtained from the so- 
lution of the corresponding homogeneous equation, 



The probability for an approaching electron in the 
continuum state k to make a transition to any of the 
bound states n of the polarization-induced image poten- 
tial is given by the prompt energy-resolved sticking coef- 
ficient,— 



prompt 



n 



(1) 



where r t = 2L/v z is the travelling time through the sur- 
face potential of width L which, in the limit L —> oo, can 
be absorbed into the transition probability from the con- 
tinuum state k to the bound state n, W nk . If the incident 
unit electron flux (we consider only a single electron im- 
pinging on the surface) is stationary and corresponds to 
an electron with Boltzmann distributed kinetic energies, 
the prompt energy-averaged sticking coefficient is given 
by^ 



prompt _ 



E 



prompt £ . 



k °e,k 



J2 k ke~P° E * 



(2) 



where /3" 1 = kgT e is the mean electron energy. 

Prompt sticking coefficients are properly weighted 
sums over state-to-state transition probabilities from 



— n n (t) = ^2 [W n n>nn'(t) - W n > n n n (t)\ - ^2w kn n n (t) 



T nn > TL n > (t) , 



(5) 



and treating the electron flux j k (t) as an externally spec- 
ified quantity . 21 ! 22 In the simplest case, which is also the 
basis of Eq. ©, jk(t) is the stationary flux correspond- 
ing to a single electron whose energy is Boltzmann dis- 
tributed over the continuum states k with a mean elec- 
tron energy k B T e , that is, j k (t) = j k ~ ke~ p " Ek . 

To solve Eq. (JSJ amounts to solving the eigenvalue 
problem for the matrix Tj 21 ' 22 For the particular case of 
an electron physisorbing at a dielectric surface this has 
been already done in I. If the transitions between bound 
states are much faster than the transitions to the contin- 
uum, so that the adsorbed electron escapes very slowly, 
one eigenvalue, — Ao, turns out to be considerably smaller 
than all the other eigenvalues — X K . The equilibrium oc- 
cupation of the bound states, n° q , is then to a very good 
approximation the right eigenvector to — Ao, which can 
be thus identified with the negative of the inverse of the 
desorption time, that is, Ao = t^ 1 . 

The kinetic sticking coefficient, which takes into ac- 
count not only the initial capture but also the subsequent 
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relaxation of the occupancy of the bound surface states, 
can be obtained as follows. — The solution of Eq. §3§ is 
split according to 



n n {t) = rf n {t) + n n {t) , 



(6) 



where 



<{t) = e- Aot f dt'e^'e^ Y, 4 0) nW lk j k (t') (7) 



is the slowly and 



ni(t) =E e ^"' / '4 K) E e f )r '^(*') ( g ) 

K#0 J -°° k,l 

the quickly varying part of n n (t). The quantities and 
e„ are, respectively, the components of the right and left 
eigenvectors of the matrix T to the eigenvalue — A K . The 
probability of the electron remaining in the surface states 
for times of the order of the desorption time is given by 
the slowly varying part only, that is, n s (t) = J2n n n(t)- 
Differentiating n s (t) with respect to time, 



the momentum transfer parallel to the surface when the 
electron absorbs or emits a surface phonon. Taking this 
effect into account leads to a recoil-corrected image po- 
tential ~ l/(z + z c ) with z c a cut-off parameter defined 
in I. Transitions between the surface states supported by 
the image potential are caused by a longitudinal acoustic 
bulk phonon perpendicular to the surface. 

The Hamiltonian from which we calculated the tran- 
sition probabilities was introduced in I where all quanti- 
ties entering the Hamiltonian are explicitly defined. It is 
given by-2^ 



where 



H = H 



static 



h, 



ph 



H 



dyn 
-ph 



qC q C q 



(12) 



(13) 



describes the electron in the recoil-corrected image poten- 
tial, which thus accounts for the coupling of the electron 
to the dipole-active surface phonon, 



H p h = 



(14) 



dt 



^^E^w-vw, (9) 



enables us, following Brenig 21 , to identify the kinetic 
energy-resolved sticking coefficient, 



kinetic 
s e,k 



(10) 



71,1 



which gives the probability for the electron being trapped 
even after the energy relaxation of the second stage of ph- 
ysisorption. In analogy to Eq. @ the energy- averaged 
kinetic sticking coefficient reads for a stationary Boltz- 
mannian electron flux 



kinctii 



s kinetic£,g-0 e .E fc 



(11) 



III. TRANSITION PROBABILITIES 

The transition probabilities W qq > , where q and q' stand 
either for fc or n, are the fundamental building blocks of 
the foregoing analysis. They have to be calculated from 
a microscopic model for the electron-surface interaction. 
The necessary steps have been described in I. 

In short, for a dielectric surface, the main source, leav- 
ing interband electronic excitations aside, which primar- 
ily affect the dielectric constant, of the attractive static 
electron-surface potential is the coupling of the electron 
to a dipole-active surface phonon . 42 > 43 Far from the sur- 
face the surface potential merges with the classical image 
potential and thus ~ 1/z. Close to the surface, however, 
it is strongly modified by the recoil energy resulting from 



describes the free dynamics of the acoustic bulk phonon 
responsible for transitions between surface states, and 

Htl h =Ytt\^z)\<l)cl,c q . (15) 

9:9' 

denotes the dynamic coupling of the electron to the bulk 
phonon. Expanding V p (u, z) with respect to the displace- 
ment field, 



E 



2/j,loqN s 



(16) 



allows us to classify the dynamic interaction according to 
the number of exchanged bulk phonons. 

As in I we use a Debye model for the bulk phonon. 
Sums over phonon momenta are thus replaced by 



E 

Q 



3iV s 
u 3 D 



dcouj 



(17) 



Measuring energies in units of the Debye energy frioo = 
^bTbi important dimensionless energy parameters char- 
acterizing (|T2"j) are 



Tiuj 



and A„„/ = 



E n — E n i 



D 



Tiuj 



D 



(18) 



where E n < is the energy of the n th bound state. We 
call the surface potential shallow if the lowest bound state 
is at most one Debye energy beneath the continuum, that 
is, t\ > — 1, one-phonon deep if the energy difference 
between the lowest two bound states does not exceed 
one Debye energy, that is, A12 > —1, two-phonon deep 
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TABLE I: Dielectric constant e s , Debye energy Hujd, en- 
ergy difference of the lowest two bound states of the recoil- 
corrected image potential AE12, and the corresponding po- 
tential depth parameter A12 for graphite, silicon dioxide, and 
gallium arsenide. 

e s fajJD Ai?i2 A12 
graphite 13.5 0.215eV 0.233eV 1.06 
Si0 2 3.8 0.041eV 0.105eV 2.59 
GaAs 13 0.030eV 0.152eV 5.13 

if the energy difference between the lowest two bound 
states is between one and two Debye energies, that is, 
— 1 > A12 > —2, and so forth. 

Because of the strong interaction between the electron 
and the dipole-active surface phonon, physisorption of an 
electron typically takes place in an at least two-phonon 
deep image potential (see Table |T| . Hence, physisorption 
of an electron controlled by a bulk acoustic phonon, as 
anticipated in (|12[) and in fact applicable to dielectric sur- 
faces, where large energy gaps block electronic relaxation 
channels due to internal electron- hole pairs and/or plas- 
mons, has to involve the exchange of many bulk phonons. 

The transition probability from an electronic state q 
to an electronic state q' is given by«iS> 

a,s' s 

x S(E S - E a , + E q - E q ,) , (19) 

where T is the on-shell T-matrix corresponding to H e ^ h 
and I3 S = (ksTs)^ 1 with T s the surface temperature; s) 
and \s') are initial and final phonon states, which are 
averaged over. 

Multi-phonon processes have two possible origins. 16 
They arise from the expansion of H^^ h with respect to 
u, 

Ht a ph =V 1 + V 2 + V 3 + 0(u 4 ) (20) 

and from the multiple action of this perturbation. Defin- 
ing the free electron-phonon resolvent, 

Go = (E-H s ^ tic -H ph + te)- 1 , (21) 

the latter is encoded in the T-matrix equation, 

T = H^ ph + H^; h G T. (22) 

Using the short-hand notation introduced in I, the one- 
phonon process, proportional to u 2 , is accounted for by 

( S / ,q'\V 1 \s,q)(s,q\V 1 *\ S , ,q') . (23) 

It leads to the standard golden rule approximation for 
the transition probability. 

Two-phonon processes are proportional to it 4 and thus 
less likely than one-phonon processes. Most of them 



TABLE II: Material parameters for the numerical results. 



Debye temperature 


T D 


2500K 


Dielectric constant 


£s 


13.5 


TO phonon mode frequency 


hiJT 


170meV 


Griincisen parameter 


7G 


1.7 


Shear modulus 


M 


5 GPa 



renormalize only the one-phonon transition probability 
and can thus be neglected in a first approximation. There 
are however two-phonon processes which induce transi- 
tions absent in the one-phonon approximation and hence 
have to be included in the calculation of the transition 
probabilities. In our short-hand notation the processes 
in question are 

(s',q'\V 2 \s,q)(s,q\V 2 *\/,q') , (24) 

( S , ,q , \V 2 \s,q)(s,q\V 1 *G*V*\ S , ,q') , (25) 

(s',q'\V 1 G V 1 \s,q}(s,q\V 2 *\s',q') , (26) 

(s',q'\V 1 G V 1 \s,q}(s,q\V 1 *G* V 1 *\s',q') . (27) 

It is shown in I how these processes can be included in the 
calculation of the transition probabilities W qq i entering 
the rate equation ([3]). Singularities appearing in some 
of the two-phonon transition rates have been regularized 
by taking a finite phonon lifetime into account (see I and 
Ref. 44 for details). 

The electronic matrix elements entering the transition 
probabilities have been also calculated in I, using bound 
and unbound wavefunctions of the recoil-corrected image 
potential. Hence, not only bound states but also contin- 
uum states belong to the static surface potential^ Our 
approach is thus on par with the distorted-wave Born 
approximation employed, for instance, by Armand and 
Manson for the calculation the sticking coefficient for 
light neutral particles.— 

IV. RESULTS 

The material parameters chosen for the numerical cal- 
culations are, unless specified otherwise, given in Ta- 
ble HIl They correspond to graphite. For some calcula- 
tions we use however the Debye temperature as a tunable 
parameter to realize different potential depths which is 
the main focus of this investigation. 

A. One-phonon deep potentials 

First, we present results for shallow and one-phonon 
deep surface potentials. In leading order, only one- 
phonon processes are involved and the one-phonon ap- 
proximation for the transition probabilities is sufficient. 
Because the electron thermalizes then very quickly with 
the surface the prompt and kinetic sticking coefficients 



5 




FIG. 1: Upper panel: Sketch of a shallow (left) and one- 
phonon deep (right) potential. The grey shaded areas show 
the energy range of sticking by one-phonon processes. Mid- 
dle panel: Energy- resolved prompt sticking coefficient as 
a function of the electron energy for a shallow potential 
(T D = 4100K) at T s = 410K (left) and for a one-phonon 
deep potential (T D = 3000K) at T s = 300K (right). Lower 
panel: Energy-averaged prompt sticking coefficient as a func- 
tion of the mean electron energy for a shallow potential 
(Tb = 4100K) at T s = 205K (left) and a one-phonon deep 
potential (To = 3000K) at T s = 150K (right). 

are almost identical. In this subsection we show there- 
fore only results for the prompt sticking coefficient. 

Figure Q] compares s^ TC " :npt for a shallow and a one- 
phonon deep potential. The sketches in the upper part 
of the figure illustrate the main difference between the 
two potentials. For a shallow potential the lowest bound 
state is less than one Debye energy below the continuum 
so that a low-lying electron from the continuum can be 
directly trapped in the lowest bound state, n = 1, by a 
one-phonon transition. In the case of a one-phonon deep 
potential, one-phonon processes can only lead to trapping 
in one of the upper bound states n > 1. 

The middle panels of Fig. Q] show the prompt energy- 
resolved sticking coefficient as a function of the energy 
of the incident electron. Apart from discontinuities the 
sticking coefficient depends linearly on the electron en- 
ergy. As explained in Sec. IIIII the one-phonon transi- 
tion probability is proportional to u 2 . From Eq. (|I6[) 
we have u 2 ~ 1/uj so that in conjunction with the Debye 
model (|17p the transition probability is proportional to oj 
which translates due to energy conservation to a propor- 
tionality to the electron energy. The phonon spectrum is 



thus reflected in the (one-phonon) energy-resolved stick- 
ing coefficient, as it is, for instance, also in the cross 
section for (one-phonon) inelastic particle-surface scat- 
tering^ 

Steep jumps in the energy-resolved sticking coefficient 
reflect the level accessibility. When the energy difference 
between the electron and a bound state exceeds the De- 
bye energy, one-phonon transitions are no longer possible 
and the electron can no longer directly reach that level. 
For a shallow potential, the first drop reflects therefore 
the accessibility of the first bound state, whereas for a 
one-phonon deep potential, where this bound state can- 
not be directly reached, the first drop reflects the acces- 
sibility of the second bound state. As energy differences 
between successive bound states of the image potential 
decrease towards the ionization threshold, that is, with 
increasing n (see upper panels of Fig. [1]), more such steps 
are found near the maximum electron energy allowing for 
trapping, which is the Debye energy. 

The contribution of the n bound state to the stick- 
ing coefficient, reflected in the height of the correspond- 
ing accessibility threshold, decreases for higher bound 
states. The reason for this lies in the electronic ma- 
trix element appearing in first order perturbation the- 
ory, (n\l/(z + z c ) 2 \k). This matrix element is smaller for 
higher bound states because higher bound states have 
less weight near the surface where the perturbation is 
strongest. Of considerable importance is hence the low- 
est bound state, which, if available, yields a particularly 
large contribution. The decreasing electronic matrix ele- 
ment also implies that neglecting all but a few, say seven, 
of the infinitely many bound states suffices for the calcu- 
lation of the sticking coefficient. 

The prompt energy-averaged sticking coefficient is 
shown in the lower panels of Fig. [T] as function of the 
mean electron energy. Due to thermal averaging the 
sticking coefficient does no longer exhibit characteris- 
tics of the phonon spectrum and level accessibility, mak- 
ing it thus more robust against changes in the phonon 
model. It does however reflect the importance of the low- 
est bound state for shallow potentials. Note also, due to 
the long-range tail of the recoil-corrected image potential 
~ 1/ (z + z c ) the energy-resolved and the energy-averaged 
electron sticking coefficients are finite for vanishing elec- 
tron energy and electron temperature, irrespective of the 
surface temperature, as it should bei 26 i 31 



B. Two-phonon deep potentials 

We now turn our attention to two-phonon deep poten- 
tials. Under the assumption that the true one-phonon 
process dominates the corrections coming from two- 
phonon processes, the latter are only taken into account 
for transitions where one-phonon processes alone would 
yield no transition probability. 

Two-phonon processes affect in a two-phonon deep po- 
tential sticking in two ways. They enable prompt trap- 
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E e [eV] 

FIG. 2: Energy-resolved prompt sticking coefficient for a 
two-phonon deep potential (T D = 2500K and T s = 500K) 
calculated with different numbers of bound states N. Full 
lines denote the one-phonon contribution, dashed lines the 
two-phonon contribution. Inset: Contribution of the second 
bound state. One-phonon contribution (red), two-phonon 
contribution (blue) broken down into the processes (V2) 2 , 
(Vi) 2 T/ 2 and (V1) 4 . 



ping from higher-lying continuum states, outside the one- 
phonon trapping range, and they control the energy re- 
laxation of the trapped electron and thus the formation of 
the quasi-stationary bound state occupancy from which 
desorption occurs. There are thus two questions to be 
answered: How significant are two-phonon processes for 
prompt sticking and how does the relaxation thereafter 
depend on the type of phonon processes available. 

To address the first question we show in Fig. [5] the 
contributions to the prompt energy-resolved sticking co- 
efficient arising from, respectively, one- and two-phonon 
processes. If available, one-phonon processes provide for 
much larger sticking coefficients than two-phonon pro- 
cesses. Figure [D also confirms that the sticking coeffi- 
cient saturates quickly with the number of bound states 
included into the calculation. 

To investigate the relative importance of the various 
two-phonon processes arising, respectively, from the ex- 
pansion of the dynamical perturbation and the T-matrix 
we plot in the inset of Fig. [5] the partial contributions 
to the prompt sticking coefficient arising from the vari- 
ous two-phonon processes which trigger transition to the 
second bound state. A two-phonon process can be simul- 
taneous, as encoded in V2, or successive, as described by 
V\GqV\. Hence, the total two-phonon transition proba- 
bility contains a contribution without virtual intermedi- 
ate states, symbolically denoted by V 2 2 [see Eq. (|24p ] and 
two contributions with virtual intermediate states, sym- 
bolically denoted by {V X ) 2 V 2 and (V1) 4 [see Eqs. (g5)), 
(|26|) . and (|27)) ]. The prompt energy-resolved sticking co- 
efficient calculated with either V 2 2 , (Vi) 2 V2, or Vy only 
is shown in the inset of Fig. [5J In accordance to what 
we found in our calculation of the desorption time of an 




0.1 0.2 0.3 0.4 0.5 

E e [eV] 

FIG. 3: Prompt and kinetic energy-resolved sticking coeffi- 
cient as a function of the electron energy for a two-phonon 
deep potential (T D = 2500K and T s = 357. 14K). Full lines 
one-phonon contribution, dashed lines two-phonon contribu- 
tion. 



image-bound electron (paper I) and to what Gumhalter 
and Siber found in their calculation of the cross section 
for inelastic particle-surface scattering^—, the direct 
two phonon process V 2 2 is dominated by the processes 
V1V2 and (Vi) 4 . 

Having clarified that two-phonon processes lead to 
a much smaller prompt sticking coefficient than one- 
phonon processes we now move on to study the effect 
of two-phonon transitions on the relaxation of the bound 
state occupancy. For a two-phonon deep potential the 
energy difference between the lowest two bound states 
exceeds one Debye energy. Hence, the relaxation of an 
electron trapped in one of the upper bound states to 
the quasi-stationary occupancy can only occur via two- 
phonon processes. Since the kinetic sticking coefficient 
gives the probability for the incident electron making 
not only a transition to a bound state but also a sub- 
sequent relaxation to the quasi-stationary occupancy of 
these states, the importance of two-phonon processes 
should be signalled by the amount the kinetic sticking 
coefficient deviates from the prompt sticking coefficient. 

Figure [3] shows that for a two-phonon deep potential 
the kinetic energy-resolved sticking coefficient is for in- 
termediate electron energies considerably smaller than 
the prompt energy-resolved sticking coefficient. This is 
due to the fact that the two-phonon transition to the 
lowest bound state, where the major part of the quasi- 
stationary occupancy resides, is very unlikely and thus 
very slow. Only for small energies, the first hump of the 
two-phonon contribution to the sticking coefficient, due 
to trapping in the lowest bound state, are prompt and 
kinetic sticking coefficients identical because no trickling 
through is needed. 

The weak coupling between the lowest two bound 
states in a two-phonon deep surface potential leads to 
a relaxation bottleneck for the electron if it is initially 
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FIG. 4: Upper panel: Inverse desorption time from the lowest 
level (dashed blue line) and the upper levels (full red line). 
Middle panel: Probability for an electron initially trapped in 
one of the upper levels of the surface potential (n = 2, 3, 4 . . .) 
either to fall to the lowest bound state (dashed blue line) or 
to desorb without ever reaching the lowest bound state (full 
red line). Lower panel: Prompt (full red line) and kinetic 
(dashed blue line) energy-resolved sticking coefficient. In all 
three panels, E e = O.leV and T s /Td = 0.2 (to keep the level 
of phonon excitation constant we set Td/T s constant). For 
Td < 2707K the surface potential is two-phonon deep, for 
2707K < T D < 4029K it is one-phonon deep, and for T D > 
4029K it is shallow. 



trapped in one of the upper states. In Figs. [4] and [5] 
we analyze the relaxation bottleneck in more detail as a 
function of the Debye temperature Tjj (to realize different 
potential depths) and the surface temperature T s . The 
upper panel shows the desorption time from the lowest 
bound state, that is, the desorption time for an electron 
capable to fall to the lowest bound state, and the des- 
orption time from the upper bound states, that is, the 
desorption time for an electron not capable to fall to the 
lowest bound state. The probability for the electron ini- 
tially trapped in the upper bound states to fall down to 
the lowest bound states and the probability to desorb to 
the continuum without ever passing through the lowest 
bound state are shown in the middle panel and the lower 
panel shows the prompt and the kinetic sticking coeffi- 
cient. 

Before we discuss Figs. H] and [SJ we say a few words 
about the way we calculated the quantities shown in the 
upper and middle panels. The desorption time from the 
lowest bound state is the desorption time for an electron 
equilibrated with the surface, the quantity we calculated 
in I, because the quasi-stationary occupancy and the 
equilibrium occupancy coincide and both reside more- 
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FIG. 5: The three panels show, as a function of the surface 
temperature, the quantities of Fig. |4]for To ~ 2500K, that is, 
graphite, and E e = 0.09eV. 



over, for the considered surface temperatures, mainly on 
the lowest level. The desorption time from the upper 
bound states we simply calculated from the rate equa- 
tion (J5J) with the lowest bound state excluded. 

The probabilities shown in the middle panels we 
obtained from the following consideration. Whether 
an electron trapped in the upper bound states passes 
through to the lowest bound state or not depends on how 
large the transition probabilities from the upper bound 
states to the lowest bound state are in comparison to the 
transition probabilities to the continuum. Hence, the sec- 
ond stage of physisorption, that is, the time evolution of 
the occupancy after the initial trapping, can be captured 
by a rate equation for the occupancy of the upper bound 
states (n — 2, 3, . . .), similar to Eq. ([5]), but with a loss 
term to both the continuum and the lowest bound state, 

-T^n n = [Wnn'n>n'(t) - W n > n n n (t)] 



- ^2 W knn n (t) - W ln n n (t) 
k 

-■D nn 'Tl n ' (£) , 



(28) 



where n and n' run over the upper image states. Solving 
Eq. (|28p with the initial condition 



nz(0) = 



Efc Wikjk 

J2l,k W lkjk 



(29) 



which is the (conditional) probability that the electron is 
trapped in the I th image state under the condition that it 
is trapped in any of the bound states, we deduce for the 
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probability for an electron trapped in one of the upper 
bound states to fall either to the lowest bound state (/ = 
1) or to desorb without falling to the lowest bound state 
(/ = c), 



p f = n f (t -^oo)=J2 Wfn^-d^ #°»l(0) 



(30) 



where, dfcf 1 and <$k are the components of the right and 
left eigenvectors of the matrix D. 

We now turn our attention to Fig.[4j It shows the effect 
of different potential depths realized by tuning the De- 
bye temperature. For a shallow potential (To > 4029K) 
desorption from the lowest level is mainly due to direct 
one-phonon transitions to the continuum; the same type 
of transition emptying the upper bound states. Hence, 
the desorption time from the lowest bound state and the 
upper bound states, respectively, differ not too much 
for shallow potentials. For one-phonon deep potentials 
(2707K < T D < 4029K), however, the cascade of two 
one-phonon processes via the second level yields much 
larger desorption times from the lowest level compared 
to the desorption time form the upper levels. For a two- 
phonon deep potential (Td < 2707K), finally, the first leg 
of the cascade, the transition to the second bound state, 
is a two-phonon transition, which increases the desorp- 
tion time compared to a one-phonon deep potential. 

The second level is the link between the upper bound 
states and the lowest bound state. The ratio of the transi- 
tion probabilities from the second bound state to the low- 
est bound state, W12, and from the second bound state 
to the continuum, W C 2 , determines if the electron trickles 
through after initial trapping or not, that is, whether it 
thermalizes with the surface or not. For a one-phonon 
deep potential both W\2 and W C 2 are due to one-phonon 
processes, in this case W12 > W C 2- For a two-phonon 
deep potential, however, the transition from the second 
to the lowest bound state is enabled by a two-phonon 
process only. In this case, and for moderate surface tem- 
peratures, W12 < W C 2, so that the electron is more likely 
to desorb before relaxing to the lowest bound state. As 
the kinetic sticking coefficient gives the probability of the 
trapped electron to relax to the quasi-stationary occu- 
pancy, the drop in the probability for reaching the low- 
est level at Td — 2707K, which is the one-phonon/two- 
phonon threshold, translates into the abrupt reduction of 
the kinetic sticking coefficient at Td = 2707K (see middle 
and lower panels of Fig. 0]). 

Figure [5] shows the quantities of Fig. @] as a function of 
the surface temperature. The Debye temperature is fixed 
to the value for graphite. At low surface temperatures 
the kinetic sticking coefficient is only slightly smaller 
than the prompt sticking coefficient, yet for high sur- 
face temperatures their difference increases significantly 
as a consequence of the inhibited thermalization. This 
can be understood as follows: The transition from the 
second to the first bound state entails the emission of 
two phonons and the transition from the second bound 
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FIG. 6: Bound state occupancy of the lowest bound state 
(dashed blue line) and the upper bound states (solid red line) 
as a function of time for an unit flux of a Boltzmannian elec- 
tron with ksTe — O.leV. The left and right panels show the 
two occupancies on two timescales. The left panel on the scale 
of the desorption time A^" 1 = r e (vertical dashed line in the 
left panels) and the right panel on the scale set by A^ 1 (ver- 
tical dashed line in the right panels). The upper two panels 
are for T s — 200K and Td = 2500K whereas the lower two 
panels show results for T s = 600K and T D = 2500K. 



state to the continuum requires only the absorption of a 
single phonon. At low enough surface temperatures it is 
however nevertheless possible that the electron drops to 
the lowest bound state because the likelihood of phonon 
emission is proportional to I + ns whilst the likelihood 
of phonon absorption is proportional to jib- Hence, for 
sufficiently low surface temperatures, W12 > W e 2, even 
when W12 entails a two-phonon and W C 2 a one-phonon 
process, so that the electron has a good chance to trickle 
through. Increasing the surface temperature benefits 
however W C 2 more than W12 so that W12 < W C 2, pro- 
hibiting the trickling through and leading to a consider- 
able reduction of the kinetic sticking coefficient at high 
surface temperatures. 

From the discussion of Figs. Q] and [5] we conclude that 
a pronounced relaxation bottleneck inhibiting thermal- 
ization can only occur for at least two-phonon deep po- 
tentials and sufficiently high surface temperatures. 

The question arises then on what timescale does the re- 
laxation bottleneck affect physisorption. To answer this 
question we analyze, respectively, the time evolution of 
the occupancy of the lowest level and the occupancy of 
the upper levels of the surface potential under the as- 
sumption that initially all bound states were empty and 
that for t > a stationary unit flux of a Boltzmannian 
electron fills the levels. Accordingly, the occupancy of the 
lowest state (n = I) and the upper states (n > 2) can be 
determined from Eq. © setting jk(t) = for t < and 
jk(t) = Jk ~ ke-P° E * for t > with J2k 3k = 1. 

The results of this calculation are shown in Fig. [5] 
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FIG. 7: Prompt (full line) and kinetic (dashed line) energy- 
averaged sticking coefficient for graphite (To = 2500K) as a 
function of the mean energy of the electron and the surface 
temperature. 



for low (upper two panels) and high (lower two panels) 
surface temperature. Clearly, for times of the order of 
the desorption time, r e — Aq" 1 , indicated by the verti- 
cal dashed line in the left panels, the upper levels are 
basically empty indicating that a thermalized electron 
desorbs; for T D = 2500K and T s = 500K the quasi- 
stationary occupancy deviates from the equilibrium oc- 
cupancy less than 3%. The upper levels are more popu- 
lated than the lower one only for very short timescales, 
set by A^ 1 , indicated by the vertical dashed line in the 
right panels. Since A^ 1 <C A^ 1 , the relaxation bottle- 
neck does not affect desorption, which still occurs from 
the equilibrium occupancy. It does thus only affect the 
kinetic sticking coefficient which is significantly smaller 
than the prompt one and actually the one to be used to 
characterize polarization-induced trapping of an electron 
at a dielectric surface. The relaxation bottleneck is ab- 
sent in neutral physisorption systems because the level 
spacing is small compared to the Debye energy. Prompt 
and kinetic sticking coefficients are thus almost identical 
as has been indeed found for neon atoms physisorbing on 
a copper substrate^ 

Figure [7] finally shows for graphite the energy-averaged 
prompt and kinetic sticking coefficients as a function of 
the mean energy of the incident electron and the surface 
temperature. As two-phonon processes contribute little 
to the initial trapping of the electron, their most im- 
portant role is to control relaxation to the lowest bound 
state. In agreement with the foregoing discussion, the ki- 
netic sticking coefficient diminishes therefore for higher 
surface temperatures whereas the prompt sticking coef- 
ficient is less sensitive to the surface temperature. From 
Fig. [7J it can be also seen that even the prompt sticking 
coefficient for graphite is only at most of the order of 
of 10 ~ 3 , the order we also found in our investigation of 
electron sticking at metallic surfaces^ It is two orders of 
magnitude smaller than the value obtained from a semi- 



classical estimate^ whose range of applicability is how- 
ever hard to grasp. We expect it, at best, to be applicable 
to very low mean electron energies, below 0.0026eV, and 
rather high electron binding energies, larger than leV^ 



V. CONCLUSIONS 

As a preparatory step towards a microscopic under- 
standing of the build-up of surface charges at dielectric 
plasma boundaries, we investigated phonon-mediated 
temporary trapping of an electron on a dielectric sur- 
face. In our simple model for the polarization-induced 
interaction of the electron and the dielectric surface, the 
adsorbed electron occupies the bound states of a recoil- 
corrected image potential. Electron energy relaxation re- 
sponsible for transitions between the image states leading 
to adsorption and eventually to desorption is due to the 
coupling to an acoustic bulk phonon. 

Dielectrics typically used as plasma boundaries are 
graphite, silicon oxide, aluminium oxide, and bismuth 
silicon oxide. They all have large energy gaps blocking 
internal electronic degrees of freedom and small Debye 
energies compared to the energy difference of at least the 
lowest two bound surface states. Electron physisorption 
at these boundaries is thus driven by multi-phonon pro- 
cesses. As in I we presented results for a two-phonon deep 
surface potential, as it is applicable to graphite, where the 
energy difference between the lowest two bound states is 
between one and two Debye energies. Classifying two- 
phonon processes by the energy difference they allow to 
bridge, we included two-phonon transition probabilities 
only for transitions not already triggered by one-phonon 
processes. Besides the Debye temperature, which we var- 
ied to realize different potential depths, the material pa- 
rameters used in the numerical calculations are the ones 
for graphite. 

Similar to physisorption of a neutral particle, stick- 
ing and desorption of an electron can be subdivided 
into three characteristic stages. At first, the electron is 
trapped in one of the upper bound states of the surface 
potential. Then the bound state occupancy relaxes to a 
quasi-stationary occupancy. Finally, over the timescale 
set by the desorption time, the electron desorbs. In or- 
der to account for both initial trapping and subsequent 
relaxation we employed a quantum-kinetic rate equation 
for the occupancy of the image states. Apart from calcu- 
lating the energy-resolved and energy-averaged prompt 
and kinetic electron sticking coefficients, which typically 
turn out to be of the order of 10 -3 , we also investigated 
the relative importance of one- and two-phonon processes 
for the two stages of the sticking process. 

The initial trapping is almost entirely due to one- 
phonon transitions from the continuum, two-phonon pro- 
cesses from higher-lying continuum states contribute very 
little. The relaxation of the bound state occupancy af- 
ter the initial trapping depends strongly on the ratio of 
the probabilities for downwards transitions to the low- 
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est state and upwards transitions to the continuum. For 
graphite, with its two-phonon deep surface potential, the 
upper bound states are linked to the lowest bound state 
only by a two-phonon process. The trapped electron has 
thus only a slim chance to drop to the lowest bound 
state, particularly at high surface temperatures, which 
favor transitions back to the continuum. The decreased 
accessibility of the lowest surface state leads to a sig- 
nificant reduction of the kinetic sticking coefficient com- 
pared to the prompt sticking coefficient. For the other 
dielectrics typically used as plasma boundaries, silicon 
dioxide, aluminium oxide, and bismuth silicon oxide, the 



surface potentials are much deeper because the Debye 
energy for these materials is very small. Hence, more 
than two phonons are required to link the upper image 
states to the lowest one, the accessibility of the lowest im- 
age state is thus even more suppressed, and the kinetic 
sticking coefficient should be accordingly small. 
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